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ABSTRACT 

We explore the dynamical and thermal evolution of the ejected neutron star crust in a Quark- 
Nova explosion. Typical explosion energies and ejected crust masses result in relativistic ejection 
with Lorentz factors of a few to a few hundred. The ejecta undergoes a rapid cooling and stretching 
resulting in break up into many small pieces (clumps) when the ejecta is only ~ 100 km from the 
explosion site. The number and size of the clumps depends on whether the breakup occurs in the 
liquid or solid phase. For these two cases, the clump number is ~ 10 3 (liquid phase) or ~ 10 7 (solid 
phase) and, at break up, are spherical (size ~ 10 4 cm; liquid phase) or needle shaped (~ 10 4 x 10 2 
cm; solid phase). 

Subject headings: Collapsar - QuarkNova - Supernova - quark star - GRBs 



1. INTRODUCTION 

In the quark-nova (QN) picture, (Ouyed et al. 2002; 
Keranen et al. 2005; hereafter ODD and KOJ re- 
spectively) the core of a neutron star, that undergoes 
the phase transition to the quark phase, shrinks in a 
spherically symmetric fashion to a stable, more compact 
strange matter configuration faster than the overlaying 
material (the neutron-rich hadronic envelope) can re- 
spond, leading to an effective core collapse. The core 
of the neutron star is a few kilometers in radius initially, 
but shrinks to 1-2 km in a collapse time of about 0.1 ms 
(|Lugones et al. 19941) . The gravitational potential en- 
ergy released (plus latent heat of phase transition) dur- 
ing this event is converted partly into internal energy and 
partly into outward propagating shock waves which im- 
part kinetic energy to the material that eventually forms 
the ejecta. 

There are three previously proposed mechanisms for 
ejection of the outer layers of the neutron star (i.e. 
crust) : (i) Unstable baryon to quark combustion leading 
to a shock-driven ejection (Horvath & Benvenuto 1988). 
More recent work, assuming realistic quark matter equa- 
tions of state, argues for strong deflagration (Drago et al. 
2007) that can expel surface material. In these models 
up to 10~ 2 M Q can be ejected. These calculations focus 
on the microphysics and not on the effect of the global 
state of the resulting quark core which collapses prior to 
complete combustion (KOJ), leading to conversion only 
of the inner core (~ 1-2 km) of the neutron star; (ii) 
Neutrino-driven explosion where the energy is deposited 
in a thin (the densest) layer at the bottom of the crust 
above a gap separating it from the collapsing core (KOJ). 
For the neutrino-driven mechanism, the core bounce was 
neglected, and neutrinos emitted from the conversion to 
strange matter transported the energy into the outer re- 
gions of the star, leading to heating and subsequent mass 
ejection. Consequently, mass ejection is limited to about 
10~ 5 M Q (corresponding to the crust mass below neu- 
tron drip density) for compact quark cores of size (1-2) 
km; (iii) Thermal fireball driven ejection which we con- 
sider for the present study. The fireball is inherent to the 
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properties of the quark star at birth. The birth temper- 
ature was found to be of the order of 10-20 MeV since 
the collapse is adiabatic rather than isothermal (ODD; 
KOJ). In this temperature regime the quark matter is in 
the superconducting Color-Flavor Locked (CFL) phase 
(Rajagopal&Wilczek 2001) where the photon emissivity 
dwarfs the neutrino emissivity (Vogt, Rapp, & Ouyed 
2004; Ouyed, Rapp, & Vogt 2005). The average pho- 
ton energy is ~3T for a CFL temperature T: since the 
plasma frequency of CFL matter is Tilo ~ 23 MeV (e.g. 
Usov 2001), the photon emissivity is highly attenuated 
as soon as the surface temperature of the star cools be- 
low T a = huip/3 ~ 7.7 MeV. To summarize, the thermal 
fireball is generated as the star cools from its birth tem- 
perature down to ~ 7.7 MeV. 

Here we focus on the case where a QN goes off in iso- 
lation, neglecting any interaction with the surroundings, 
and the crust is ejected by a thermal fireball (i.e. case 
(iii) above). The remainder of this paper is presented 
as follows: In Sect. 2 we discuss the energetics of the 
QN and the resulting fireball and ensuing crust ejection. 
Sect. 3 deals with the hydrodynamical and thermal evo- 
lution of the ejecta. We conclude in Sect. 4. 

2. FIREBALL GENERATION AND CRUST EJECTION 

The QN ejecta, which is the left-over crust of the par- 
ent neutron star, is initially in the shape of a shell and is 
imparted with energy from the QN explosion. This en- 
ergy released during the QN explosion (-Eqn,53; in units 
of 10 53 erg) is a combination of baryon to quark con- 
version energy and gravitational energy release due to 
contraction. Also, the rapid contraction creates a gap be- 
tween the surface of the collapsed quark matter core and 
the inner edge of the remaining hadronic matter. Before 
entering the CFL phase at T ~ 10-20 MeV the QS cools 
mainly by neutrino emission; as shown in Keranen et al. 
(2005) a significant fraction of neutrinos escape through 
the crust. Once the QS core enters the CFL phase, the 
energy is released in photons (see Vogt et al. 2004; Ouyed 
et al. 2005). Then at that point thermalization leads to 
the creation of a (e + e _ ) fireball with a temperature of 
10-20 MeV. Thus the initial state for crust ejection is a 
QS core surrounded by a fireball in the gap (~ 1 km) in 
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turn surrounded by the initially stationary neutron star 
crust. 

The general picture is that the fireball from the QN 
acts as a piston at the base of the crust; i.e. the QS fire- 
ball expands approximately adiabatically while pushing 
the overlaying crust, and cooling fairly rapidly. Using 
the relativistic force equation, a simple estimate for the 
timescale to accelerate the crust to Lorentz factor of a few 
is ~ 0.1 to 1 millisecond. The force is due to the fireball 
pressure (F = PA where P = (a/3)T 4 , T ~ 10-20 MeV 
and A is the shell's area). The acceleration timescale is 
only a factor of a few smaller than the free-fall timescale. 
As the crust moves outward, the fireball cools approx- 
imately adiabatically, with the temperature T propor- 
tional to V^ 1 / 3 with V the volume between the quark 
star surface and the inner edge of the crust. Initially 
the temperature decreases slowly with R(t) then later it 
decreases as l/R(t) where R(t) is the radial position of 
the crust. Furthermore, the final kinetic energy of the 
crust for a Lorentz factor of 10 is about 2 x 10 50 erg for 
a 10~ 4 M© crust. Thus the energy needed to eject the 
crust is less than 1% of fireball energy. 

The energy input to the crust is mainly kinetic and 
only a small fraction is thermal since there is no strong 
shock in this scenario. During the acceleration, the crust 
can respond on timescale AR/c s where AR is the crust 
thickness (about 1 km) and c s is the speed of sound 
in the crust. The majority of the crust is degenerate 
with c s = c/a/3- The response timescale is ~ 5 x 10~ 6 
seconds, about (l/100)th of the acceleration timescale. 
Thus no strong shock will be generated except in the 
outermost non-degenerate (i.e. low sound speed) layers 
(the mass contained in the non-degenrate layers is small 
< 10~ 6 Mq). The bulk of the crust is thus accelerated 
smoothly by the pressure from the fireball. In compari- 
son, a normal SN has the energy of core collapse released 
into mostly non-degenerate matter. Thus a strong shock 
is generated. In the QN case, most of the crust is highly 
degenerate with high sound speed, thus it can be accel- 
erated without a strong shock. 

One can estimate the Lorentz factor 1 of the QN ejecta 

by 
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where 770.1 is the efficiency of energy transfer from the 
QN to the ejecta's kinetic energy in units of 0.1, and the 
ejecta mass, m C j OCta ,-4 is given in units of 10 -4 M Q . Mass 
ejection in the QN scenario depends on details of energy 
transfer to the crust. Best estimates of ejected mass from 
existing calculations are 1O -5 Mq-1O~ 2 M0. We adopt a 
rounded, fiducial, ejected mass of 10 _4 M Q for the re- 
mainder of this paper. Heavier mass ejection would lead 
to mildly or non-relativistic ejecta. 

3. EVOLUTION OF THE QUARK-NOVA EJECTA 

For ejecta masses higher than 10 the correspond- 

ing initial density and size (i.e. of the neutron star 

1 For simplicity we assume that the entire ejected crust can be 
represented using a single Lorentz factor, density and temperature. 
However, a more realistic situation might involve ejecta with a dis- 
tribution of these variables. This can only be determined with 
calculation of the ejection process which requires complex simula- 
tions. 



crust) just before ejection are po > 10 10 g cm™ 3 and 
Ar > 0.025r respectively (e.g. Datta et al. 1995), 
where ro is the radius of the parent neutron star before 
the QN explosion. Also, the ejecta's birth temperature, 
To, is estimated to be of the order of 10 MeV. Using 
general volume expansion (V = Vo(r/ro) Qv ), mass con- 
servation implies that the density of the ejecta will be 
p = po(r/ro)~ av . We define ay by a power law depen- 
dence of V on r with different ay describing different 
physical situations, e.g. constant volume, adiabatic, etc.. 
At this density and temperature the ejecta is relativistic 
degenerate, so that the Fermi energy evolves with radius 
as ep oc p 1 / 3 oc (^) _ " v ^ 3 (see Appendix IA1 for details). 

The total internal energy, U, in the ejecta at birth can 
be estimated to be Uq ~ 5 x 10 47 erg TO c jccta,-4 /°o io- 
Here, the ejecta mass is given in units of 10~ 4 M Q , its 
density in units of 10 10 g cm -3 , and fi e = 2 is the mean 
molecular weight dominated by the degenerate electrons. 
The corresponding heat equation describing the thermal 
evolution of the ejecta can then be written as 



dU __ 
~dt ~ L 
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dt 
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where all quantities are expressed in the ejecta's frame. 
Since the ejecta is expanding relativistically, in its frame 
we can approximate dr — cdt with r in this case being 
the distance to the star from the ejecta. This allows the 
equation above to be recast into, 
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3.1. Relativistic expansion of a degenerate ejecta 

We solve for the heat equation (Eq. [3]) for the cases 
of adiabatic and isothermal expansion of the QN ejecta 
in the shape of a spherical shell. The expansion is rela- 
tivistic and assumed to occur at at constant speed. The 
analysis is in the ejecta's frame. In the adiabatic case, 
L 7 = which leads to U/U = (r/r )~ 2Qv/3 . Using the 
expression for tp (with e^o — £f{po)) given in Appendix 
A we arrive at 
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which shows that the temperature goes to zero in finite 
time corresponding to a radius 
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Clearly the radius is too close to the origin which means 
that in the adiabatic case the ejecta would cool almost 
immediately and would continue to expand at zero tem- 
perature, with U — (3/5)Nep. In reality the ejecta will 
solidify when it cools below the melting temperature of 
iron (see §2.2 below). 
In the isothermal case, T = To , the heat equation leads 

to 

U ={ r L 2av/3 _ A 
U r 
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( r o,6 T o 4 io/ m -4Pio 3 ); r o,6 is r Q in units 
of 10 6 cm and T i is T in units of 10 MeV. The internal 
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energy U goes to zero at a finite radius 
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which corresponds to r ~ 1.15ro for a v ~ 9/4. The 
corresponding ejecta density at this radius is ~ 10 9 gm 
cm -3 . For T ~ 0.1 MeV we hnd r ~ 17r for a v ~ 9/4. 
with a corresponding ejecta density ~ 10 8 gm cm -3 . 

In summary, applying the heat equation to the rela- 
tivistically expanding ejecta shows an almost instanta- 
neous loss of the internal energy. This mainly due to the 
PdV work done by the ejecta as it expands or due to its 
rapidly increase in area which leads to efficient cooling. 

3.2. Liquid to solid transition 

In the early stages the QN ejecta, with a tempera- 
ture in the tens of MeV, resembles a hot molten plasma. 
Because it is degenerate it cools extremely rapidly (see 
Appendix A). As the ejecta moves outwards it expands 
and cools undergoing a liquid to solid transformation - 
we will refer to the corresponding radius as the solidifica- 
tion radius r s . The state of ions (in this case iron nuclei) 
can be conveniently specified by the Coulomb plasma pa- 
rameter (e.g. Potekhin et al. 1999), 
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where Z and A are the ion charge number and the atomic 
weight of iron. Since here 2 > 1 prior to solidification, 
the shell material constitutes a strongly coupled liquid. 
Solidification occurs for 3 > S m = 172 which implies a 

1 /3 

solidification temperature of T s ~ 9.5p 8 keV (see also 
de Blasio 1995). For example, using equation above and 
for a v ~ 9/4, an ejecta born at ~ 10 MeV and po — 10 10 - 
10 11 g cm -3 , will cool to below T s ~ 10 keV and solidify 
(see §2.1) when it reaches a radius of ~ 3-10ro. The so- 
lidification radius is then r s ~ 30-100 km, at which point 
the ejecta's density is p s > 10 8 g cm -3 (see Appendices 
13. ip . The corresponding ejecta density at this radius is 
~ 10 8 gm cm~ 3 . 

Previous studies of decompressed neutron star crust 
show that the material heats up as it expands because 
the matter fuses into heavier elements and releases en- 
ergy, which could prevent crystallization. The relevant 
calculations for expanding quark-nova ejecta has been 
performed in detail in Jaikumar et al. (2007). That work 
considers r-process in quark-nova ejecta with and with- 
out /3-decay heating. In particular section 3 in that paper 
discusses the decompression of the ejecta. As can be seen 
from Figures 5 and 6 in Jaikumar et al. (2007) the ejecta 
does not reheat above 10 keV (see peak in temperature 
at time of a few milliseconds) . The temperature prior to 
/3-decay reheating is low enough for crystallization and 
break-up to occur. However, /3-decay reheating which 
peaks at a few milliseconds after ejection (see figures 5 
and 6 in Jaikumar et al. 2007) might re-melt the crystal- 
lized ejecta. Even so, the melting temperature and the 
temperature of the re-heated ejecta may be similar: this 
would lead to crystallization at some later time. Finally, 
we note that those calculations did not include radia- 
tive losses so they provide an upper limit on the ejecta 
reheating temperature, thus increasing the likelihood of 
crystallization. 



TABLE 1 

Clump/chunk properties at breakup for Liquid (L) and 
Solid (S; see Appendix cases 



10" 3 190 



1CT 1 1.9 



if> Ar c (cm) g c (rad) m c (gm) jV c 

190 2 x 1(T 5 1.7 x 10 19 2.5 x 10 
x~10 4 2 x 10~ 3 1.7 x 10 23 2.5 x 10 



2.5 x 10 3 



3.3. Ejecta Breakup 

Regardless of when crystallization occurs, the relativis- 
tic expansion causes rapid breakup into small chunks be- 
cause of the inability of causal communication laterally 
in the shell. Specifically, as the spherical shell expands 
radially outwards relativistically, it also stretches later- 
ally relativistically 2 . Since the lateral expansion of the 
matter in the shell is limited by the speed of sound, c s , 
adjacent patches of the shell separate leading to breakup 
into small clumps. 

The breakup depends on the existence of surface ten- 
sion in the expanding shell. For liquid or solid iron (our 
case), inter-ionic forces (mediated by the electrons) pro- 
vide the tension. Such a tension does not exist in a gas. 
E.g., for SN ejecta where the shock accelerates the gas 
beyond the sound speed the lack of surface tension means 
the particles simply expand away from each other. An 
analogy is that of a bursting balloon versus expanding 
dust. For the innermost SN ejecta the density may be 
high enough that it behaves like a liquid rather than a 
gas. In that case, it is likely the sound speed (which in- 
creases inward) is higher than the expansion speed in the 
inner parts (which decreases inward; w tx r). If the SN 
ejecta is still a liquid at the radius where the expansion 
speed exceeds the speed of sound then it will be subject 
to this lateral breakup mechanism. 

The size of the clumps depends on whether the breakup 
occurs in the liquid or solid phase. We discuss the solid 
phase case first; the liquid case differs in having a much 
larger breakup strain ip. As the ejecta solidifies at r s , 
a strain rapidly builds up inside as the ejecta continues 
to expand relativistically. Defining the breaking strain of 
solid iron to be -0 ~ 10 -3 (e.g. Halliday & Resnick §13.6) 
then the ejecta break-up occurs at roughly a radius of 
Tb ~ (1 + i>)r s (see Appendix ICj) . The typical chunk size 
at birth is c s At where At = (n, — r B )/c, or, 



186 cm Cb^b. 

1/6 



II Ml 



(9) 



where we defined (b = ip-3P\, g ■ Here ij) and the breakup 
radius are given in units of 10~ 3 and 100 km respectively, 
while the density at solidification radius is in units of 10 s 
g cm -3 . 

Table 1 lists the properties of the clumps/chunks for 
the two cases of breakup occurring in early liquid or later 
solid phase. If the shell breaks into clumps during the 
liquid phase, there are two possibilities: 

• (i) the chunks will breakup into smaller pieces when 
they solidify. Because of the degeneracy of the 
ejecta, the cooling is so rapid that even if initial 
breakup occurs in the liquid phase it will be imme- 
diately followed by another breakup governed by 
the smaller breaking strain ip of the solid. Since the 
breakup radius is hardly changed the solid chunk 

2 The lateral separation velocity for 2 points separated by A9 is 
cA8 (see Appendix [Ct. 
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size should be independent of whether clumping 
first took place during the liquid phase or not. 

• (ii) the ambient temperature remains high enough 
to keep the clumps liquid as they expand. We note 
that as the density drops due to expansion, the 
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(see Appendix [Cl for details). Furthermore, their length 
can be found, from equation (|B2|) . to be Ar c j ccta .scp. ~ 
5 x 10 6 cm mV^ ta _ 4 . Hence, we expect the chunk's 
length is roughly a hundred times its width for case (i) 
or equal to its width for case (ii). In the observer's frame 
the ratio of length to width will be contracted by a factor 
of 1/IY For case (i), the chunks thus resemble what we 
will refer to as "iron needles" , reminiscent of the subjet 
model of Toma et al. (2005) in the context of GRBs. Fi- 
nally, as the chunks expand radially outwards to a radius 

r 2 

r ut > ?"scp, we can associate a filling factor / out = - J r E -. 



crystallization temperature also drops (T m oc r 



-2/3 

since p cx r~ 2 ). There is a delicate competition 
between the dropping ambient temperature and the 
dropping crystallization temperature. 

In summary, if during the radial expansion, the am- 
bient temperature drops below T s , then the liquid will 
crystallize. We adopt case (i) as our fiducial case. How- 
ever most of our results are the same for both cases and 
we note whenever case (ii) gives a different outcome. In 
what follows we use the term chunks to refer to either 
the small solid iron pieces or the larger liquid clumps. 

3.4. Expansion of broken ejecta 

Beyond the break-up radius, r^, the chunks remain 
in contact with each other within the relativistically ex- 
panding ejecta as a whole. This is because the pieces 
expand in volume, filling up the space between them, 
and causing the density of each piece to continuously de- 
crease, until they reach the zero pressure iron density 
(pFc ~ 10 g cm~ 3 ), at which point they stop expanding. 
If we define a filling factor to be, /b = 1 at rb, then 
it will remain unity until some separation radius r scp .. 
This radius of separation can be found from equation 
(|B3[) . With p = py c and our fiducial values this radius is 



' scp. 
?"0 



1300 
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which implies that the chunks remain closely packed until 
the ejecta as a whole reaches r. 



scp. 



10 cm. This is an 
upper limit for the separation radius since we neglected 
the effects of the ions (see Appendix B). 

An estimation of the cross-sectional area extended 
by each chunk just before separation is given by 



4. SUMMARY AND CONCLUSION 

In this paper we investigate the thermal and dynamic 
evolution of a relativistically expanding iron-rich shell 
from a QN explosion. The QN produces a photon fire- 
ball which acts as piston to eject and accelerate the crust 
of the parent neutron star to mildly relativistic speeds. 
We find that the shell rapidly cools and breaks up into 
numerous (10 3 to 10 7 ) chunks because of the rapid lat- 
eral expansion. Breakup may occur while the material is 
in liquid or solid phase: if liquid the chunks are nearly 
spherical; if solid they are needle shaped moving parallel 
to their long axis. 

Although the presented model is based on physical 
arguments, most of these are in reality more compli- 
cated and so would require more detailed studies and the 
help of numerical simulations. For example, the process 
of clumping, crystallization, and breakup of the ejecta, 
would require better knowledge of the ambient conditions 
surrounding the ejecta. The astrophysical implications 
will be discussed elsewhere. 
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APPENDIX 



FERMI GAS 

The Fermi energy for a relativistic and a non-relativistic gas is ep 

a 2 /3 2/3 

1.219 x 10 erg p 10 ~ 76.2 MeV p 10 respectively (e.g. Shapiro&Teukolsky 1983; p24). The internal energy, u 



per particle is 



1.413 x 10- 5 erg p^ 3 



8.8 MeV p\£ and 
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with a corresponding pressure P = (2/3)(u/v) where v is the volume per particle. In the relativistic degenerate regime, 



io 15 ^ 473 ; 



P = K r p ' with k t = 1.244 x 
The corresponding sound speed is c s , r 



in the non-relativistic regime P = n nl .p 5 / 3 with K r 



2.57 x 10V /6 



and 



2.26 x 



10V /3 



transition from relativistic to non-relativistic degeneracy occurs at density ptr — 2 x 10 6 g cm 
from Fermi-Dirac to Boltzmann statistics occurs at the degeneracy temperature kT ~ ep- 



cm s 

3 



9.91 x 10 12 p~ 5/3 . 
x , respectively. The 
Finally, the transition 



EVOLUTION OF EJECTA THICKNESS AND DENSITY 

We approximate the pressure in the ejecta by degenerate electron gas pressure neglecting the contribution from the 
non-relativistic ions. The ejecta thickness increases in time at the speed of sound d(Ar) — c s dt with dt — dr/c in 
the ejecta's frame. It is straightforward to show that the ejecta's density quickly drops below the transition density 
p tr ~ 2 x 10 6 g cm -3 so that most of the ejecta thickness expansion occurs during the non-relativistic degenerate 
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phase. That is, d(Ar) ~ 2.26 x 10 6 p 1 / 3 dt with p = m c j ccta / (inr 2 Ar) . Combining these equation gives 



where B — n\^ cta _ 4 t^q / Ar$^ with the ejecta's thickness at birth Aro,4 given in units of 10 4 cm, corresponding to 
po ~ 10 10 g cm -3 (e.g. Datta et al. 1995). For r » r , 

Ar ~ 2.6 x 10 4 cm mV 4 cta ,_ 4 r 1/4 . (B2) 



Mass conservation then gives (for r >> ro), 



p 10~ 2 fr Y /4 
Po 



EJECTA BREAK-UP 

A spherical shell expanding radially at a speed of light gives 2 points separated by A9 a lateral separation velocity 
cA9. This causes a strain to rapidly build up leading to breakup. Here we discuss the case that the shell has solidified 
prior to break up (see §2.1 for the liquid case). A typical breaking strain of iron is ip ~ 10~ 3 (e.g. Halliday & Resnick 
§13.6). The matter in the ejecta expands at the speed of sound (c s ) yielding a the break-up radius of, 



rb = r s 1 + ip 



r B (l + V0 , (CI) 



P c 1/4 _ 3 /4 
p Cs m cjocta,-4 r 

where r s is the solidification radius and j3c is the velocity of the ejecta. 
The typical chunk size at birth is c s At where At = (rb — r s )/c, or, 

Ar c ~ 186 cm tpspl^bAOO = 186 cm Cbn^ioo > ( C2 ) 

where we have defined (b = ipspl/g • H ere ip and the break-off radius are given in units of 10~ 3 and 100 km respectively, 
while the density at solidification radius is in units of 10 8 g cm -3 . 

The angle subtended (measured from the QS) by a broken chunk of the ejecta at birth (i.e. at the break-up radius 
r b ) is 

9 C = Ar c /r b - 2 x 10' 5 rad Cb • (C3) 

The corresponding solid angle is then, 

n c = tt9 2 - 5 x 10" 10 sr Cb • (C4) 
The corresponding mass of the broken-off chunk is then roughly, 

m c = ^m cjcct a ~ 1-7 x 10 19 gm Cb^cjccta.-4 ■ (C5) 

The ejecta break-up only depends on the stress in the two dimensions due to the ejecta's lateral expansion being 
independent of its thickness. Thus, the chunk mass is linearly dependent on the ejecta mass as a more massive ejecta 
would be thicker. The total number of chunks can be estimated as, 



3 2 T 2 9 2 
JV c = ^~2.5x 10 7 — |M , (C6) 



where the collimation angle in the ejecta frame is 9b, c = T#b- 
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